Mem. S.A.It. Vol. xx, 1 

© SAIt 2002 Memorie deiia 



Structural and scaling properties of galaxy 
clusters: probing the physics of structure 

formation 

M. Arnaud, G. W. Pratt* and E. Pointecouteau 

CEA/DSM/DAPNIA/Service d'Astrophysique, CE Saclay, L'Orme des Merisiers, 91191 
Gif sur Yvette, France 



Abstract. We present XMM-Newton studies of the total mass, gas density, temperature and 
entropy profiles in nearby hot and cool clusters, together with follow-up observations of 
distant clusters from the SHARC Survey. The observed structural and scaling properties are 
compared with the predictions of the self-similar model of cluster formation. These data 
indicate that clusters do form a self-similar population down to low mass and up to high 
redshift, and give support to the standard picture of structure formation for the dark matter 
component. However, deviations from the standard scaling laws confirm that the specific 
physics of the gas component is still insufficiently understood. 
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1. Introduction 

1.1. The self-similar model of cluster 
formation 

In the standard hierarchical formation sce- 
nario, clusters of galaxies are forming in the 
recent cosmological epoch (z ~ 2 to the 
present time). In theory. X-ray cluster for- 
mation and evolution is simple and driven 
by the collisionless gravitational collapse 
of the main dark matter (DM) component, 
whereas the hot gaseous intra-Cluster Medium 
(ICM) 'follows' the gravitational potential 
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of the DM. Analytical models of the DM 
collapse ( Bertschinaer "198?; Cava liere et al. I 
fl^ 9). sup ported by numerical simulations 
llEvrard & G ioia 2002 for a review), then pre- 
dict that galaxy clusters constitute a self- 
similar population. 

In such a model, purely based on grav- 
itation, the internal shape of clusters is uni- 
versal, independent of mass and redshift. 
Furthermore, the virialised part of a cluster, 
present at a given redshift z, corresponds to a 
fixed density contrast, 6^ ~ 200, as compared 
to the critical density of the Universe at that 
redshift, pdz)'- 



M2OO 47T 

—i— = -r^c Pc(z) 



R 



(1) 



200 



where M200 and R200 are the virial mass and ra- 
dius. pc(z) = /i^(z)3i/2/(87rG), where /z(z) is the 
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Fig. 1. Universal profiles obtained from numerical 
simulations of the formation and evolution of galaxy 
clusters (from Navarro et al. 1995). The radius is 
scaled to the virial radius, corresponding to an over- 
density of 200 as compared to the mean density of 
the Universe in which the cluster is embedded. The 
Dark matter and gas densities are normalized to the 
mean density. 



Hubble constant normalized to its local value. 
h{z) depends on the cosmological density pa- 
rameter Qm and the cosmological constant A: 
h\z) = + z)' + A. 

The ICM is assumed to be an ideal gas 
evolving in the gravitational potential of the 
dark matter. The gas mass fraction /gas is thus 
constant: 
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and self similarity applies to both the DM com- 
ponent and the ICM. The virial theorem then 
gives: 



Gjump M20Q 
2 ^200 



(3) 



where T is the mean temperature and /Sj is 
a normalisation constant depending on cluster 
(universal) internal structure. 

From the basic equations Eq.^ Eq.|2]and 
Eq. |3] it is easy to see that the cluster popu- 
lation is a two parameter population. Simple 
scaling laws relate each physical property, Q, 
to the cluster X-ray temperature and redshift, 
in the form Q oc A(z)T". The evolution factor. 



A(z), is due to the evolution of the mean DM 
(and thus gas) density, which varies as the crit- 
ical density of the Universe: 

<^ Pdm o^pciz) « h^(z) 



Pi 



(4) 



For instance, the total and gas mass scales as 
Mgas M200 °^ h^\z)T^^^, the virial radius as 
/?20() °^ h^^{z)T^^^ and the X-ray luminosity, 
assuming bremsstrahlung emission, as Lx oc 

1.2. Insights from X-ray observations 

How the statistical properties of the observed 
cluster population compare with these theoret- 
ical predictions is a key issue in modern cos- 
mology. Such comparisons provide unique in- 
formation on the physics that governs the for- 
mation and evolution of large scale structure, 
both for the DM and the baryonic components. 

We can first examine relations between 
global properties such as luminosity, gas mass, 
total mass, size and temperature. This is a 
powerful method, and has been m uch used 
with data from previous satellites llArnaud I 
12002! for a review and references below). 
Although strong coiTelations are observed be- 
tween various quantities, it was rapidly rec- 
ognized that clusters deviate from the classic 
self-similar model. The most remarkable de- 
viation is the local Lx-T relation, which is 
steeper than expected (e.g. Arnaud & Evrardl 
il999HMarkevitch 119981) . This was the first in- 
dication that the gas physics should be looked 
at more closely and that non- gravitational pro- 
cesse s could play a role (e.g. lEvrard & Henry I 
l1991h . 

An unambiguous picture can be obtained 
only by studying the internal structure of clus- 
ters. For instance a steepening of the Lx-T re- 
lation could be due to a systematic increase of 
/gas with T (a simple modification of the scal- 
ing laws) and/or a variation of cluster shape 
with T (a break of self-similarity). By look- 
ing at radially averaged profiles of interesting 
quantities (temperature, density, entropy, inte- 
grated mass, etc) two considerations can be ad- 
dressed at the same time. The first is to see 
whether the profiles agree in shape; the second 
is to investigate the sealing of the profiles. The 



Arnaud, Pratt & Pointecouteau: Structural and scaling properties of galaxy clusters 



3 



existence of a universal underl ying dark mat- 
ter pr ofile in the simulations (iNavarro et al. I 
Il997t hereafter NFW) spurred searches to de- 
tect the observable signatures of such univer- 
sality in the gas density and temperature pro- 
files. This was indeed found to be the case 
when ROSAT and ASCA data of hot clusters 
were examined ( e.g. Markevitc h et al. lll998l: 
iNeuman & Am aiid 1999). ROSAT/ASCA ob- 
servations also gave the first indication that the 
self-similarity does h old up to high redshift 
jArnaud et al. l2002al) . 

The present generation of X-ray satellites 
represents a giant step forward in terms of res- 
olution and sensitivity. In marked contrast to 
studies with the previous generation of satel- 
lites, we are subject to far less uncertainty on 
the temperature profiles (and thus on the total 
mass and entropy profiles) due to PSF effects. 
Furthermore, the sensitivity of XMM-Newton 
now allows us to probe the scaling and struc- 
tural properties further away from the clus- 
ter centre, down to lower mass clusters and at 
higher precision at high redshifts. 

2. Dark matter in nearby clusters 

2.1. A universal dark matter profile 

The observed shape of the DM distribution is 
an important test of our understanding of the 
DM collapse. We derived dark matter profiles 
from the XMM-Newton observations of three 
nearby clusters of galaxies. Two of these are 
hot and massive: A1413 (z - 0.143; mean 
temperature YT - 6.5 keV; iPratt & Arnaud I 
[2002) and A4 78 (z ^ 0,0881; kr = 6.7 keV; 
[Pointecoutea u et al. 1 120031) . The third cluster 
is the cool system A1983 (z = 0.044), with 
kr = 2.3 keV (Pratt & Arnaud 2003). The to- 
tal mass profiles were obtained from the tem- 
perature and gas density profiles under the hy- 
pothesis of hydrostatic equilibrium and spher- 
ical symmetry. Data processing details (back- 
ground subtraction, vignetting and PSF cor- 
rection, deprojection) can be found in the pa- 
pers. The total mass profiles are well fitted us- 
ing a dark matter model as described by NFW. 
In the case of A478 the radial coverage was 
wide enough (especially in the centre) to dis- 
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Fig. 2. Integrated total mass of A1413 (black di- 
amonds; Pratt & Arnaud 2002), A478 (red trian- 
gles; 'Pointecouteau et al. 2003) and A1983 (blue 
squares; Pratt & Arnaud 2003) scaled to their virial 
mass versus the scaled radius (XMM-Newton data). 
The best fit NFW model profile is overplotted in the 
same color code for each cluster. 



tinguish between the NFW and the more cen- 
tral ly peaked profile fou nd in the simulations 
by iMoore et al. I (Il999h . Such a profile was 
found to be inconsistent with our data, as well 
as profiles without a cusp, like the King profile. 

The three profiles were scaled according 
to their virial mass, M200, corresponding to a 
density contrast of 200, and their virial radius 
{R200), as shown in Fig. |2 The three scaled 
profiles are in good agreement. From these 
three XMM-Newton observations, we are able 
to describe the shape of a universal dark mat- 
ter profile in nearby clusters in the radial range 
~ 0.01/?2{)o to 0.77? 200- Our result suggests that 
this profile shape is independent of the cluster 
temperature (or mass) as predicted in the hier- 
archical model of structure formation. 

2.2. TheM-T relation 

A further test of the scaling properties of 
nearby clusters is provided by the M-T re- 
lation. No definitive answer could be derived 
from A5CA/7?(9&4r observations, various stud- 
ies providing diff'erent scalings in both normal- 
isation and slope. For instance, it was unclear 
whether the mass sc ales as T^^^ as expected 
("Horne r et al. 1119991) . or if this was true only 
in the high mass regime {kT > 3 - 4 keV, 
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Fig. 3. M - T relation at /?2500 as seen by 
XMM-Newton (red data points and solid line) 
and Chandra (black data points and solid line, 
[Allen et al. 2002). The blue dotted line is the p re- 
diction from numerical simulations ilEvrard et al. 



see also lEttori et al. Ij 2002t). wi th a steepen- 
ing at lo wer m ass (|N evalainen et al. | [2 00(1 
IXu et al. I '2001', 'Finogueno v et al. I l200lh .~or 

even if the slope was higher all over the mass 
range (Sande rson et al. 2003). This is likely to 
reflect both the large uncertainties in the tem- 
perature profiles (see below), and various ex- 
trapolation biases, especially for low mass sys- 
tems. 

We computed the total mass of our three 
clusters at a radius corresponding to a density 
contrast of 2500 (M2500). A power law fit was 
performed in the M - T plane, leading us to 
an observational M - T relation as seen by 
XMM-Newton: M2500 This rela- 

tion can be compared with the relation derived 
at_the same density contrast by Allen et al. 
J2OO2I) from Chandra observations of six hot 
and massive clusters: M2500 <^ 7'i-52±o.36 gQj-jj 
results are in excellent agreement (Fig. [^ 
and are consistent with the standard scaling. 
The precision on the slope obtained with the 
XMM-Newton data from a sample of only three 
clusters is remarkable. This is both due to the 
precision of the individual mass estimates and 
the unique leverage provided by XMM-Newton 
at the low mass end. 

However, comparison with the predic- 



terms of normalisation. This d iscrepancy ap- 
pears to be present at all radii ( Pratt & Arnaud I 
l2002h . Knowing that the normalisation of the 
M -T relation depends on both the gas density 
and dark matter profiles, this result strongly 
suggests that, while the dark matter collapse 
seems to be correctly understood, we are still 
unable to model correctly the gas profile in 
terms of shape and distribution. In the next sec- 
tion we will discuss the gas properties. 

3. Some insights into cluster gas 
pliysics 

3.1. Emission measure profiles 

It is first useful to consider the emission 
measure along the line of sight EM{r) - 



"°° Kg dl, which is easily derived from the 
X-ray surface brightness profile and is directly 
linked to the gas content and density distri- 
bution. The scaled EM profiles of hot clus- 
ters were considered by Neumann & Arnaud 
(1999; 2001; see also Vikhlinin et al. 1999) 
and were found to be similar in shape outside 
the cooling core region. Recall that in the stan- 
dard self-similar framework, the mean gas den- 
sity does not depend on the temperature (see 
Eq. [5. The emission measure thus scales as 
7^200^ i-e EM oc r"^. These authors found that 
the scatter was considerably reduced if a much 
steeper EM — T relation, EM oc r'-^^, was 
instead used. This explains the steepening of 
the Lx-T relation and translates into Mgas k 
r'^'*, consistent with th e observed steepen ing 



also Vikhlinin et al. 
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Castillo-Morales & Schindler 


2003) as 


com- 



tion from numerical simulations Evrard e t al. I 
(kl996.) exhibits a significant discrepancy in 



pared to the standard Mgas o^T scaling. 

In Fig.[4[we s how the EM profile of A I983 
(kr = 2.13 keV; Pra tt & Arnaud]l2003h com- 
pared to the mean profile for hotter clusters 
(kr > 3.5 keV; the envelope indicates typi- 
cal dispersion about the mean). In the left-hand 
panel, the profiles are scaled using the stan- 
dard self-similar scaling, EM oc r"'^. In the 
right hand panel, the empirically-derived scal- 
ing of EM oc 7" 1-38 used. The use of the non- 
standard scaling not only reduces the scatter 
in the mean profile for hotter clusters, it also 
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Fig. 4. The XMM-Newton scaled emission measure (EM) profile of A 1983 (data points) compared to the 
mean scaled EM profile of the ensemble of hot clusters (kr > 3.5 keV) observed by ROSAT and studied by 
iNeuman & Amaud1i200ll) . The latter is indicated by the dashed line; dotted lines show the dispersion. On 
the left, the EM profiles are scaled using the self-similar scaling EM oc r"-^ (correspon ding to Mp„s oc r'-^). 
On the right, it is scaled using the empirical EM oc r'^^ (M^.^^ oc t'--'^*) scaling of lNeuman & ArnaudI 
<200lh. 



brings the EM profile of A1983 into the enve- 
lope of values for hotter clusters. 

This result is interesting because it implies 
a universal gas density profile shape down to 
quite low temperature — the EM profile of 
A1983 is not flatter than that of hotter clusters. 
The result also shows that the observed EM — 
T relation is indeed steeper than expected, con- 
sistent with the observed steepening of the 
Mcas — T relation. 



3.2. Temperature profiles 

The exact form of cluster temperature pro- 
files is a topic which has generated a huge 
amount of debate in the literature. Depending 
on the study, cluster profiles a re steeply de- 
cHning ( Markevitch et al. 119981) [A5CA], have 
a isothermal c ore beyond which there is a 
steep decline llDe Grandi & Molendi' "200?) 
[Beppo-SAX]), or are flat to the limits of detec- 
tion (White 2000 rA.9CA1: llrwin & Bregmanl 
[2000 [Beppo-SAX]). The differing conclusions 
reflect the complexity properly of treating the 
energy-dependent PSF of these satellites. The 
PSF of XMM-Newton is far less of a problem. 



being far smaller, and much less energy depen- 
dent. 

In Fig. Is] we show the scaled tempera- 
ture profiles of a haphazard sample of clus- 
ters which have been worked on by the Saclay 
group. Of these clusters, two display merger 
signatures (Coma, A2163), while the rest seem 
relatively relaxed. The average temperatures 
range between 2.1 and 14 keV, while the red- 
shifts lie between z = 0.02 and z = 0.6. Despite 
the differing properties, the scaled temperature 
profiles are remarkably similar. What disper- 
sion there is, is connected to the clear, well re- 
solved drop in the central regions, the signature 
of a cooling core. The profiles are essentially 
isothermal (within ~ +10%) up to ~ 0.5/? 200, 
after which there may be a decline, although 
we are presently limited by the small number 
of data points in this region. The great advan- 
tage of XMM-Newton is, of course, that we ac- 
tually have temperature measurements this far 
out from the cluster centre. 

3.3. Entropy profiles 

The entropy, defined, as is now customary, as 
S = kT In^J^, is a fundamental characteristic of 
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Fig. 5. Compilation of scaled kr profiles observed 
by XMM-Newton. The temperature is scaled to 
the mean temperature outside the cooling flow re- 
gion. The data on Coma, A1795, A2163, RXJ1120, 
A1835, A1413, and A1983 are from Arnaud et al. 
(2001a); Ai-naudetal. ( 2001b); Pratt et al. (2001); 

[Arnaud etal. (2002b); Maierowicz et al. (2002), 
Pratt & Arnaud ( 2002), Pratt & Arnaud ( 2003), re- 
spectively. The grey area corresponds to a variation 
of ±10% around unity. 



the ICM, because it is a probe of t he thermo- 
dynam ic history of the gas (e.g I Voit et al"! 
l2Q03h . In the standard self-similar picture, the 
entropy, at any scaled radius, should scale sim- 
_ply as S PC h(z )*'^T. The pioneering work s of 
[Ponman et al. I lll999h andlLlovd-Davies e t al. I 
( t200Qi) have largely been superseded by the re- 
cent work offPonman et al. ( 2003), where they 
study the entropy profiles and scaling of 66 
nearby systems observed by AS'CA and ROSAT. 
These authors conclude that, while there is 
considerable dispersion particularly at the low 
mass end, the entropy measured at 0.1/?20() in 
fact scales as 5 oc with cooler systems 

having systematically more entropy than ex- 
pected. This 'entropy excess' has been the sub- 
ject of a large amount of debate in the lit- 
erature. Proposed explanations have included 
heating before or after colla pse (from SNs or 
AGNs), radiative cooling . . . dVoit et al. Il2002t 
[Borgani et al. 2002 and references therein). 

In Fig. |6l we show the XMM-Newton en- 
tropy profiles of four clusters. Three (A1983, 
A 1991 and A27I7) are cool clusters (kr ~ 
2.0-2.5 keV; Pratt & Arnaud, in prep.); the last 
is the massive cluster A1413 {kT - 6.5 keV; 
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Fig. 6. Scaled XMM-Newton entropy profiles using 

5 cx (1 +z)-2r''-«5. 



iPratt & Arna ud "2002^. The profiles have been 
scaled using the empirically determined S oz 
jQ.65 j-gjafjon. These are the first high quality 
profiles which range from 0.01-0.57?2()o, and 
show the remarkable similarity of shape' and 
the excellent normalisation of the profiles with 
the adopted S — T relation. The external slope 
of the entropy profiles is roughly consistent 
with the S oc r' ' behaviour expected from 
shock heating. 

These profiles allow us to rule out simple 
preheating models for the entropy excess. As 
shown in Tozzi & Norman (2001) and other 
works, preheating the gas leads inescapably 
to lower mass systems having large isentropic 
cores. It can be seen that the scaled core sizes 
of all of these systems is similar, and that there 
is no evidence for a large isentropic core in the 
cooler systems. 



4. Evolution of clusters 

The evolution of cluster properties contains 
key information to i) help disentangle the re- 
spective role of gravitational and various non- 
gravitational processes, and ii) to fully test the 
self-similar model of the DM collapse. 



' We note that recent results from Chandra sug- 
:est far more dispersion in groups at about 1 keV 
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Fig. 7. Correl ation between bolornetric lumi nosity and ternperat ure (fini = 0.3, A = 0.7). Blue lines: local 
relations from I Amaud & Evrard ] il999l) andlMark evitchl \\99s ) for comparison. Left panel: Data points 
for the medium z (green points) and high z (red points) SHARC sample observed with XMM-Newton. All 
the data points stand significantly above the local relation, indicating evolution. Middle panel: Luminosities 
scaled by /j(z)"', removing the evolution expected in the standard self-similar model. There is a slight 
indication that the evolution is actually stronger . Right panel : Lum inosities have been scaled by (1 + 
z)~^^^, according to the best fit evolution found bv lLumb et al. I i200 3V Black points: Chandra data from 
IVikhlininetaiT<2002ft . Note the excellent agreement between Chandra and XMM-Newton. 



The standard self-similar model makes 
definitive predictions on the evolution of clus- 
ter properties. While they should keep the same 
internal structure as nearby clusters, more dis- 
tant clusters should be denser, smaller and 
more luminous (see Sec. 1). ROSAT/ASCA 
observations gave the first indication that the 
self-similarity does hold at z > 0. The 
universal EM profile appears to extend to 
z ~ 0.8, with a redshift scaling consistent 
with the expectation f or a ACDM cosmology 
(lArnaud et al. 1 120023). These authors found 
significant evolution in the normalisation of 
the Lx-T relation, consistent with the self- 
similar model (as have other studies consid- 
ering this cosmolo gy; iReichart et al. 1 1 19991 
lNovickietani2002 h.' However these observa- 
tions were highly biased towards massive sys- 
tems, mostly clusters discovered by the EMSS, 
and their statistical quality was poor. 

Several X-ray samples of distant clus- 
ters have been assembled using ROSAT ob- 
servations. They are much larger and cover a 
much wider mass range than the EMSS. With 
XMM-Newton it is now possible to make de- 
tailed studies of these clusters. 

We present here results from a follow- 
up of the clusters detected in the Bright and 
South SHARC survey (iRomer et al. . .2000: 



iBurke et al. 112003 '). Eight high z (0.45 < z < 
0.62) clusters w ere observed in Guaranteed 
Time (Lumb et al. 2003, Arnaud et al. 2002b) 
and six medium redshift (0.3 < z < 0.4) clus- 
ters in Open Time. A combined analysis of the 
whole sample is in pro gress (Arnaud et al. , in 
prep.; see also lMaierowicz ■2003i) . 

These data definitively confirm a positive 
evolution in the normalisation of the Lx-T 
relation (Lumbetal. 2003, Fig. In the 
left panel of Fig. we show the luminosi- 
ties for the combined SHARC sample versus 
temperature. All the data points stand signif- 
icantly above the local relation. 'Lumb et al. I 
(|2203) found that the normalisation of the 
Lx-T relation for the high z sample scales 
as (1 -H f)^-^'^^o.26^ consistent with the evo- 
lution found by Vikhlinin et al. (2002) from 
Chandra data ((1 -i- z)'-^='"-^). The right panel 
illustrates this evolution: when each luminos- 
ity is scaled by (1 -H z)"^^, the data points of 
the whole XMM-Newton and Chandra samples 
become consistent with the local relation. Note 
that this best fit evolution is stronger than the 
h(z) evolution predicted by the standard self- 
similar model (see middle panel), although the 
effect is probably not very significant (work in 
progress). 
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Fig. 8. The scaled EM profiles of medium and high 
redshift SHARC clusters (green and red lines) ob- 
served with XMM-Newton, compared to the mean 
scaled profile of nearby hot clusters (full line, the 
grey area is the scatter around the mean, see Fig.|4j- 

We show in Fig. |8] the corresponding EM 
profiles, scaled according to the standard evo- 
lution: EMsc = hizT^T-^-^^EM (the empir- 
ical slope of the EM-T relation derived by 
iNeuman & Arnaudi (1200 Ih . is assumed not to 
evolve). The profiles are traced up to similar, 
or even larger, radii than the mean profile of 
nearby clusters observed with ROSAT. There 
is good agreement between individual profiles 
and the local mean profile (taking into ac- 
count the local dispersion and errors); the self- 
similarity of form up to high z is confirmed. 
However there is a hint of a systematic discrep- 
ancy: most high z profiles lie above the above 
local curve (the effect is stronger for the high- 
est z sample). This is again an indication for 
stronger evolution, as found from the Lx-T re- 
lation. 

5. Conclusions 

The picture that emerges from the 
XMM-Newton observations presented here is 
that clusters do form a self-similar population. 

The self-similarity of shape for the gas 
component (density and T profiles beyond the 
cooling core, entropy profiles) extends down 



to surprisingly low temperature and up to high 
redshift. We begin to have evidence that the 
dark matter profile is universal and presents a 
central cusp, as predicted by numerical simu- 
lations. Scaling laws do relate various physi- 
cal properties with z and T, with a slope of the 
M-T relation as expected. This supports the 
validity of the current approach for the mod- 
eling of the dark matter collapse. 

However, the gas scaling properties are 
significantly different from expected in the 
simple model. The steepeer EM-T and shal- 
lower S-T relations for nearby clusters is con- 
firmed at all radii and down to low tempera- 
ture, and are consistent with the steepening of 
the Mg-is-T and Lx-T relations. As all these 
quantities are related to the overall gas content, 
Pgas, this suggests that the fundamental devi- 
ation from the simplest model is an increase 
of p„2j. with T (or mass), or, equivalently, an 
increase of the gas mass fraction with T (if 
M oc r'-^). There is also some indication of 
stronger evolution than expected. Furthermore, 
current adiabatic numerical simulations do not 
predict the coiTect normalisation of the local 
M-T relation, suggesting that the gas profile 
shape is not correctly reproduced. 

These results confirm that the physics gov- 
erning the baryonic component is still far from 
understood. Recent modeling, including gas 
cooling and/or non gravitational heating, has 
shown that these processes can affect the scal- 
ing relations. However, the self-similarity of 
shape is a strong new constraint, and pure pre- 
heating models can already be ruled out. 

These results must still be taken with cau- 
tion. In spite of the quality of the data, they are 
based on a small data set of clusters, and should 
be considered as 'test cases'. Observations of 
large samples, with unbiased mass and redshift 
sampling, is essential i) to confirm the self- 
similarity of form, ii) to firmly establish the 
slope and normalisation of the M-T relation 
and the intrinsic scatter in cluster properties, 
and iii) fully to assess cluster evolution, in- 
cluding possible evolution of the scaling laws 
slopes. 

Acknowledgements. MA thanks D. Lumb and the 
XMM-Q project col laboration for provi ding the Sx 
profiles published in lLumb et al. I i2003h . 
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